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Abstract. MAGIC has been performing optically triggered Target of Opportunity 
observations of flaring blazars since the beginning of its scientific operations. The alerts of 
flaring blazars originate from Tuorla Blazar Monitoring Programme, which started the optical 
monitoring of candidtate TeV blazars in 2002 and has now collected up to eight years of data on 
more than 60 blazars. These ToO observations have resulted in the discovery of five new VHE 
7-ray emitting blazars (S5 07f6+7I4, IES 1011+496, Mrk 180, ON 325 and B3 2247+381). 
In addition part of the discovery of BL Lac and the discovery of 3C 279 was made during a 
high optical state. In this contribution we present a detailed analysis of the optical light curves 
which are then compared to MAGIC observations of the same sources. We aim to answer the 
question: "Is there a connection between optical and VHE 7-ray high states in blazars or have 
we just been lucky?" 

Keywords: blazars, observations, multi-wavelength 



1. Introduction 

Blazars are the most extreme case of Active Galactic Nuclei (AGN) in which the relativistic jet 
originating from close to the central black hole is oriented close to our line of sight. In these 
objects the dominant radiation component comes from the jet and is non-thermal in origin. 
Their Spectral Energy Distribution (SED) is characterised by two broad peaks of which the 
lower energy one is believed to originate from synchrotron emission of the charged particles in 
the jet. The higher energy one is most commonly explained by inverse Compton scattering of 
either the synchrotron (Synchrotron Self Compton - SSC, see e.g. [1]) or external (External 
Compton - EC, [2], [3]) seed photons by the electrons and positrons in the jet, but hadronic 
models (see e.g. [4], [5]) are also considered. Blazars are divided into fiat spectrum radio quasars 
and BL Lacertae objects. The latter are further subdivided into low-, intermediate- and high- 
frequency peaked BL Lac (LBL, IBL and LBL, respectively) objects according to the frequency 



of the first peak which typically ranges from near-infrared in LBLs to UV or soft X-rays in 
HBLs. 

Blazars show a variable flux in all wave bands from radio to Very High Energy (VHE, here 
defined as E>100 GeV) 7-rays on time scales ranging from months to a few minutes. The 
relationship between the variability at different wave bands appears to be rather complicated. 
In the simplest case both SSC and EC mechanisms predict a connection between optical and 
7-ray flares. However, there could be time lags between the flares in the different wave bands 
due to several mechanisms (see e.g. [6], [7] for further discussion) with the optical leading the 
one in 7-rays or vice versa. Also, in case the synchrotron peak is located far from the optical like 
in the case of HBLs the variability in the optical band could be very small and the synchrotron 
flares should be best visible in other wavelengths, usually X-rays. This seems to be the case for 
the best studied TeV blazars Mrk 421 [8] and Mrk 501 [9]. In principle, hadronic models can 
also explain correlated variability and cannot be excluded as the emission mechanisms of 7-rays. 

MAGIC is performing ToO observations of sources in a high flux state in other wave bands 
(e.g. optical and/or X-rays and/or high energy 7-rays). Alerts of optical high states are sent by 
the Tuorla Blazar Monitoring Program 1 [10]. Observations triggered by optical high states have 
so far led to the discovery of 5 new sources at VHE 7-rays. These are, in order of discovery: 
Mrk 180 [11], 1ES 1011+496 [12], S5 0716+714 [13], B3 2247+381 [14] and ON 325 [15]. In 
addition the discovery of 3C 279 was made and part of the data in the discovery of BL Lac were 
collected during an optical high state. Also, several 7-ray flares of known VHE emitters (e.g. 
the 2007 flare of 3C 279 [16] and the 2011 flare of 1ES 0806+524 [17]) have been detected during 
optically triggered observations. In this contribution we present an analysis of the optical light 
curves and compare them to simultaneous MAGIC observations. We discuss the implications of 
the results to the optical to VHE 7-ray connection. 

2. Analysis of the optical light curves 

Photometric observations are performed on the 35 cm Celestron telescope attached to the body of 
the 60 cm KVA telescope, which is located on the Canary island of La Palma, Spain and operated 
remotely from Tuorla Observatory, and the Tuorla 1 m telescope. Observations are done in the 
R-band and the magnitudes are measured using differential photometry with calibration stars 
in the same CCD frame as the object (see [19] for a more detailed overview). The light curves 
are available on the project web page 1 and are updated on a daily basis. 

The monitoring started in 2002 and the original list consisted of 24 candidate TeV blazars 
chosen from [18] (with 5 > 20° for them to be observable from Tuorla). Over time the number 
of monitored sources has grown to over 60 including many candidates detected at high energy 
(HE) 7-rays by the Fermi telescope. The goal of the monitoring is to study the optical variability 
of these sources and to determine their baseline flux. We have also studied the host galaxy's 
contribution to most objects' overall brightness [19]. This way we can get a direct measurement 
of the AGN's core brightness by subtracting the host galaxy contribution. For this study we 
chose 31 of the best sampled light curves including all 24 from the original sample. 

Historically a source was said to be flaring in the optical if its core flux increased by > 50% 
compared to its quiescent level. For determining the quiescent level a general rule that a source 
would spend ~20% of the time at or below this level was followed. Both the quiescent states 
and the flares were identified by visual inspection of the light curves. One of the purposes of this 
study was to calculate exact values for the quiescent states according to some common principles 
and to automate the identification of flares. As quiescent state we used the weighted average 
calculated from all the data points except for those belonging to 'flares'. For the purposes of 
these calculations 'flares' were all data points that were typically more than three (the exact 



1 http://users.utu.fi/kani/lm 



value varied slightly from source to source) standard deviations above the average. As an initial 
value for the calculations we used the historical value of quiescent state i.e. the flux of which 
20% of the whole data set was of lower level. In most cases the quiescent level ended up around 
the value where we would have put it by eye and also usually close to the 20% value used 
previously. The relative increase in core flux required for a brightening to be constituted as a 
flare was derived for each source independently depending on the amplitude of their variability 
- the larger the amplitude the higher increase in flux was required: FlareCond ~ Q state * ^ ax 

^ vain 

where FlareCond is the flare condition, Qstate the quiescent state and F max and F m i n are the 
highest and lowest measured core flux respectively. A total of 132 flares were counted since the 
beginning of the monitoring of which 108 occurred after MAGIC started regular observations. 
These numbers are slightly higher than the number of flares identified by eye from the light 
curves, because: a) if there is a long flare but with a gap in the data, software identifies it as two 
distinct flares and b) some sources (e.g. OJ 287) show "flaring states" where the flux can drop 
below the flaring condition briefly but immediately increase again which the software counts as 
separate flares. We also calculated the rise and decay times of the flares and flare durations 
where possible (refer to Fig. 1 for explanation of these terms). A summary of all the flares of 
all the sources can be seen in Table 1. 

3. Optical to VHE 7-ray comparison 

MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) is a system of two imaging 
atmospheric Cherenkov telescopes (IACTs) located on the Canary island of La Palma, Spain at 
~2200 m asl. The first telescope started observations in 2004 and since November 2009 the two 
telescopes have been observing in stereoscopic mode. Both telescopes have 17 m mirror dishes - 
the largest of any IACTs - allowing them to reach lower energy threshold than other instruments 
of their type. In the stereoscopic mode MAGIC reaches an energy threshold of ~50 GeV (~25 
GeV with a special trigger), sensitivity of ~0.8% of Crab Nebula flux (5 a in 50 h >250 GeV), 
an energy resolution of ~15-20% and an angular resolution of 0.07 degrees [20] making it the 
best instrument of its kind, particularly for energies below 200 GeV. 

Because MAGIC is performing ToO observations of blazars that are in a high state in the 
optical there is simultaneous or quasi-simultaneous optical and VHE data for quite many flares. 
Of the 108 optical flares that occurred after the beginning of Cycle 1 in MAGIC (June 2005) for 
56 there is some simultaneous MAGIC data but for only 49 of those there is more than 3 hours 
of good quality data. This is mostly constrained by the weather and moon conditions as well 
as conflicts in the schedule. On 27 occasions the source was detected at VHE (on 3 occasions 
in a low state) - a total of 14 different sources were seen (note that for Mrk 501 the VHE flare 
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Figure 1. An example of the parts of the 
flare: peak of the flare is the point in the flare 
with the highest flux; beginning and end of the 
flare are the points with the lowest flux prior or 
after exceeding the flare condition respectively; 
rise time is the time from the beginning of the 
flare to the peak and decay time the time from 
the peak to the end. 



Source 


# flares 1 


Qstate 5 


below qstate 5 


avg rise t 1 


avg decay t 5 


avg dur* 


t betw fT 


duty cycle 8 


1ES 0033+595 


3 


0.18 


15.7 


11.06 


18.98 


23.71 


73.6 


0.07 


1ES 0120+340 


1 


0.11 


20.5 


10.96 


2.99 


13.95 


- 


0.02 


RGB 0136+391 


2 


1.77 


18.1 


3.96 


55.95 


57.93 


- 


0.09 


RGB 0214+517 





0.42 


21.2 


- 


- 


- 


- 





3C SEA 


10 


5.46 


19.4 


20.96 


16.74 


36.02 


36.14 


0.22 


1ES 0647+250 


11 


0.93 


10.3 


15.3 


24.23 


35.71 


42.84 


0.38 


S5 0716+714 


9 


12.83 


23.6 


11.49 


10.3 


20.51 


45.56 


0.16 


1ES 0806+524 


3 


1.72 


31.3 


33.41 


6.02 


37.42 


19.73 


0.09 


OJ 287 


20 


2.4 


12.2 


12.08 


14.72 


24.72 


47.61 


0.33 


1ES 1011+496 


6 


2.05 


17.7 


21.8 


20.65 


42.45 


20.98 


0.22 


1ES 1028+511 


3 


0.43 


22 3 


29.49 


26.35 


46.01 


23.97 


0.13 


Mkn 421 


10 


11.41 


12.5 


16.17 


13.19 


29.36 


42.73 


0.21 


RGB 1117+202 


1 


0.58 


23.9 


45.97 


10.01 


55.99 


- 


0.06 


Mkn 180 


2 


1.13 


21 


20.65 


3.69 


24.34 


- 


0.04 


RGB113E+E76 





0.2 


29.3 


- 


- 


- 


- 





OH 325 


7 


2.27 


12.8 


14.7 


11.4 


24 


42.12 


0.22 


1ES 1218+304 


3 


0.83 


28.9 


29.14 


11.32 


30.75 


19.68 


0.13 


ON 231 


4 


2.81 


8.2 


34.31 


22.3 


42.46 


50.25 


0.35 


3C 279 


3 


1.45 


12.8 


19.95 


32.21 


45.51 


- 


0.16 


RGB 1417+257 





0.38 


20.9 


- 


- 


- 


- 





1ES 1426+428 


1 


0.41 


24.9 


14.99 


4.02 


19 




0.03 


1ES 1544+820 





0.59 


18.6 













PG 1553+113 


4 


9.31 


19.7 


34.12 


30.12 


64.24 




0.33 


Mkn 501 


4 


3.78 


19.5 


13.76 


8.18 


21 94 


55.07 


0.06 


H 1722+119 


7 


1.73 


16 


24.67 


13.96 


36.64 


40.98 


0.49 


OT 546 


1 


0.82 


21.7 


29.96 


31.81 


61.78 




0.07 


1ES 1959+650 


3 


2.98 


39 


32.84 


32.78 


65.62 




0.13 


PKS 2155-304 


5 


14.81 


18.2 


15.38 


15.14 


26.97 




0.38 


BL Lac 


6 


3.56 


19.2 


10.45 


10.34 


20.79 


69.87 


0.09 


B3 2247+381 


3 


0.47 


18.6 


40.7 


12.32 


49.93 


36.44 


0.26 


1ES 2344+514 





0.53 


26 














Table 1. Results of the light curve analysis source- by-source. 

1 Absolute value of flares identified since the beginning of the programme. 

2 Quiescent state flux in mJy. 

3 Percentage of data points below quiescent state. 

4 Average rise time in days. 

5 Average decay time in days. 

6 Average flare duration in days. 

7 Average time between two flares in the same observation period (the ~6 months when the source is 
visible). 

8 Estimate of the duty cycle (fraction of time the source spends above flaring condition). 

Note that for the calculations of the averages only values >0 were used. For the calculation of the duty 

cycles all data gaps longer than three weeks were excluded. 

occurred ~1 month after the optical), 7 of them for the first time by MAGIC. In 19 cases the 
observations resulted in upper limit only and for 8 flares information was not available or the 
analysis is ongoing. The remaining two observations resulted in the detection of VHE 7-rays 
from the 3C 66A/B region [21]. The results of the comparison of the optical and VHE data are 
summarised in Table 2. 

There are only two sources, Mrk 421 and Mrk 501, that can be reliably detected at low 
flux states in reasonable time scales and that have good enough VHE coverage that allows 
us to study the correlation between optical and VHE 7-ray states. The long-term monitoring 
of these sources by MAGIC is summarised in [22]. From the Mrk 421 data set they did not 
find any correlation between the optical and VHE wave bands whereas Donnarumma et al [23] 



source 


source type 


redshift 


VHE det?' 


# of optical flares' 


MAGIC observations' 


supports connection? 


1ES 0033+595 


HBL 


0.086* 


NO 


3 


1: -3 h. no detection 


inconclusive 


1ES 0120+340 


HBL 


0.272 


NO 


1 


no observations 


no info 


RGB 0136+391 


HBL 


>0.3 


NO 


2 


no observations 


no info 


3C GGA 


IBL 


0.444* 


YES* 


6 


5: 1 non-det, 2 det (1 of which A/B), 2 analysis ongoing 


inconclusive 


1ES 0647+250 


HBL 


0.203* 


NO 


9 


3 analysis ongoing 


inconclusive 


S5 0716+714 


LBL 


031* 


YE ST 


7 


4: 2 det, 2 analysis ongoing 


yes 


1ES 0006+524 


HBL 


0.138 


YE ST 


3 (more like 1 long) 


1: detection 


yes 


OJ 287 


LBL 


0.306 


NO 


18 


2: no detection 


no 


1ES 1011+496 


HBL 


0.212 


YES* 


6 


3: all detections 


yes 


1ES 1028+511 


HBL 


0.239 


NO 


3 


3: no detections. 2 obs > 3 h 


no 


Mkn 421 


HBL 


0.031 


YES- 


7 


3: twice in high state, once in low 


inconclusive 


RGB 1117+202 


HBL 


0.139 


NO 


1 


no observations 


no info 


m Kn 1 D U 


HBL 


0.045 


Tto 


2 


2: 1 detection, only 1 obs > 3h 


inconclusive 


ON 325 


LBL 


0.130 


YESr 


6 


3: 2 obs > 3 h, 2 non-det 1 det 


inconclusive 


1ES 1218+304 


HBL 


0.182 


YES 


2 


no observations 


no info 


ON 231 


IBL 


0.102 


YES 


4 


2: both obs <3 h. no detections 


inconclusive 


3C279 


FSRQ 


0.536 


YES- 


2 


2: both detections 


yes 


1ES 1426+42E 


HBL 


0.129 


YES 


1 


1: poor quality data no detection 


inconclusive 


PG 1553+113 


HBL 


>0.3 


YESr 


4 


2: both showed increased VHE flux 


yes 


Mkn 501 


HBL 


0.034 


YES 1 


1 


1: qamma-ray flare detected -1 month after optical 


inconclusive 


H 1722+119 


HBL 


0.018 


NO 


6 


4: no detections. 3 obs <= 4 h 


no 


OT546 


HBL 


0.055 


NO 


1 


1: no detection, strong moonlight 


inconclusive 


1ES 1959+650 


HBL 


0.048 


YESr 


3 


3: 1 obs < 2 h - no detection, twice in low state 


no 


PKS 2155-304 


HBL 


0.116 


YESr 


4 


1: high state 


yes 


BL Lac 


LBL 


0.069 


YESr 


3 


3: 1 det 2 non-det 


inconclusive 


B3 2247+381 


HBL 


0.119 


YESr 


3 (more like 1 lonq) 


1: detection 


yes 



Table 2. Summary of optical flares detected since MAGIC started observations and their comparison 
to contemporaneous VHE data. The redshifts marked with an asterisk are tentative values. 
1 Flag whether the source has been detected at VHE 7-rays by MAGIC or any other IACT. The ones 
seen by MAGIC are marked with a cross. 

2 Number of optical flares since the beginning of Cycle 1 in MAGIC. 

3 The number and short summary of the contemporaneous MAGIC observations. 

saw a hint of correlation between the two during the 2008 flare. For both Mrk 421 and 501 
extensive data sets have been collected (see e.g. [24], [25]) allowing to study the correlation 
between optical and VHE 7-rays in unparalleled detail. However, both of these sources are 
so-called extreme blazars [26] so they are not really representative of the whole blazar or the 
HBL classes. The long-term monitoring of 1ES 1959+650 by MAGIC has always detected the 
source in a low state [22] whereas in the optical the source has shown three flares since the 
beginning of the monitoring. Although difficult to observe from the Northern hemisphere due 
to large zenith angles, PKS 2155-304 can still be reliably detected by MAGIC during strong 
outbursts. One such was reported by H.E.S.S in July 2006 (coinciding with an optical flare) 
triggering MAGIC observations and leading to a strong detection of the source (see [27] and 
[28]). Also, H.E.S.S. saw a strong correlation between the optical and VHE bands during the 
multi-wavelength observations of PKS 2155-304 in August-September 2008 during a low VHE 
state of the source (see [29]) but this correlation does not always seem to be there. The results 
of five years of observations of PG 1553+113 by MAGIC reported in [30] suggest a correlation 
between the optical and VHE wave bands for this source. Only for a handful of the other sources 
we have MAGIC observations in both high and low flux states (e.g. S5 0716+714, 1ES 0806+524 
and 3C 279) and in most cases there seems to be a connection between activities in the two 
wavebands but it is difficult to prove with the low statistics at VHE. 

4. Discussion and conclusions 

Optically triggered ToO observations have resulted in the discovery of five new VHE 7-ray 
emitting blazars by MAGIC so far. Additionally several other discoveries and 7-ray flares have 
coincided with optical high states. If the synchrotron emission and the IC emission arise from the 
same emission region, a connection between optical and VHE 7-ray states would be expected. 



However, in case the synchrotron peak is located far from the optical waveband (which can be 
the case for HBLs) such connection would be difficult if not impossible to observe due to low 
amplitude of variability in the optical. And in case of different emission regions or mechanisms 
there might be no direct connection. It seems that for most sources the optical and VHE 7-ray 
states are somewhat connected. However there is at least one example (1ES 1959+650) where 
this does not seem to be the case which indicates that at least in some sources there are different 
emission mechanisms at work or that the two emissions originate from different regions. 

In order to study the problem more in detail correlation studies between the optical and VHE 
7-ray states are required. This would require more complete VHE light curves with reliable de- 
tections of the sources also at lower flux states. Unfortunately with the current generation of 
IACTs these kinds of studies are only feasible for a few sources and would require a significant 
fraction of the available observation time to be dedicated to monitoring. Otherwise we have to 
rely on future instruments to shed more light into this matter. 
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